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Abstract. 

In this letter, the integrability aspects of a generalized Fisher type equation 
with modified diffusion in (1+1) and (2+1) dimensions are studied by carrying 
out a singularity structure and symmetry analysis. It is shown that the Painleve 
property exists only for a special choice of the parameter (m = 2). A Backlund 
transformation is shown to give rise to the linearizing transformation to the linear 
heat equation for this case (m = 2). A Lie symmetry analysis also picks out the 
same case (m = 2) as the only system among this class as having nontrivial infinite 
dimensional Lie algebra of symmetries and that the similarity variables and 
similarity reductions lead in a natural way to the linearizing transformation and 
physically important classes of solutions (including known ones in the literature), 
thereby giving a group theoretical understanding of the system. For nonintegrable 
cases in (2+1) dimensions, associated Lie symmetries and similarity reductions 
are indicated. 
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The Fisher type reaction-diffusion equation with quadratic nonhnearity and 
modified diffusion of the form 

Ut-Au —(yuf-u(l-u) = 0, (1) 

1 — u 

where u{t, x, y) is some kinetic variable, A and V are gradient and Laplacian operators 
respectively and the subscript denotes partial differentiation with respect to time, is an 
important physical system appearing in many areas of physics and biology [1-4]. The 
case TO = is the classical Fisher equtaion and to = 2 corresponds to real systems such 
as bacterial colony growth [5]. It has been known for some time that the to = 2 case 
can be mapped onto the linear heat equation through a specific transformation [6]. 
Recently Brazhnik and Tyson [7] have discussed five interesting classes of travelling 
wave solutions and static structures in the 2-dimensional version of the m = 2 case 
of eq.(l). In this letter, we wish to point out that a Painleve singularity structure 
analysis [8] picks out the to = 2 case for both the (H-l) and (2-1-1) dimensions as the 
only system for which the partial differential equation (1) is free from movable critical 
singular manifolds satisfying the Painleve (P-) property. More interestingly, we point 
out that the Backlund transformation deduced from the Laurent expansion gives rise 
to the linearizing transformation in a natural way. Similarly, a Lie symmetry analysis 
singles out the to = 2 case in eq.(l) as the only system posessing a nontrivial infinite- 
dimensional Lie algebra of symmetries and that the similarity variables and similarity 
reductions give rise to the linearizing transformation and several physically interesting 
solutions, including the travelling wave solutions, static structures and so on known 
in the literature, in an automatic way. For the m ^ 2 cases, one obtains restricted 
classes of invariant solutions, including propagating pulses and fronts of special type. 

To start with, we consider the (H-l) dimensional case of eq.(l). Locally expanding 
in the neighbourhood of the non-characteristic singular manifold ^(a;, t) =0, <px,<Pt 
in the form of a Laurent series [8] 

oo 

« = 5^w,-^-+^ (2) 

j=o 

one finds that the possible values of the power of the leading order term are 

(i) P--2, 

^=1^' 
(in) p = 0. 

One can easily check that for all these three leading order cases, only for the value 
TO = 2 the solution is free from movable critical singular manifolds. One finds that in 
the case p = —1, the leading order coefficient uq is an arbitrary function in addition to 
the arbitrary singular manifold 4>{x, t). In the case of the other two leading orders, for 
TO = 2 one obtains only one arbitrary function without the introduction of movable 
singular manifolds and so they can be considered as corresponding to special solutions. 
Thus eq.(l) in the one-dimensional case is foimd to satisfy the P-property for m = 2 
and is expected to be integrable. For all other choices of to, except for certain special 
cases of the constrained singular manifold (j)t — (j)xx = 0, the solution exhibits the 
presence of movable critical singular manifold, and so the system is of non-Painleve 
type and so nonintegrable. Extending this analysis to (2-|-l) dimensional case also one 
obtains essentially the same conclusion (except for certain cases of the constrained 
manifold (pt - (pxx - <Pyy = 0). 
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Now in the Laurent expansion (2) if we cut off the series at "constant" level term, 
that is j = —p for the leading order p = 1/(1 — m) = — 1, m = 2, so that 

u=^+u^ (3) 

and demand that if wi is a solution of (1) for the case m = 2 then u is also a 
solution, one essentially has an auto-Backlund transformation. Here uq and (p satisfy 
a set of coupled partial differential equations (pdcs) arising from eq.(l) on using the 
transformation (3). Starting from the trivial solution, ui = of eq.(l), one can check 
that the equations for uq and (p in eq.(3) are consistent for the choice uq = cp, giving 
rise to the new solution u = 1 which is indeed an exact solution of (1). Now taking 
■Ui = 1 as the new seed solution, one can check from equations satisfied by uq and (p 
that 

uo = -1, (pt-(P^^- ^+1 = 0. (4) 
Defining now (p = 1 + x, one obtains the linear heat equation 

Xt-Xxx-X = 0. (5) 
Thus the transformation 

(6) 



1 + X 

where x satisfies the linear heat equation (5), is the linearizing transformation for 
eq.(l) in (1+1) dimensions for the choice m = 2. which is similar to the one considered 
in reference [6] . Here we have given an interpretation for the transformation in terms 
of the Backlund transformation. Similar analysis holds good in the case of (2+1) 
dimensions of cq.(l) also and the same transformation (6) linearizes eq.(l) for m = 2 as 
well, where x satisfies the two dimensional linear heat equation Xt ~ Xxx — Xyy ~ X = 0- 
Now let us consider the invariance of eq.(l) under one-parameter continuous Lie 
group of transformations. First we consider the (1+1) dimensional case, 

Ut - Uxx - ^ _ ^ "x -u + u =0. (7) 

Considering the infinitesimal transformation 

X — > X = X + s^{t, X, u), t — > T = t + eT{t, X, u), 

u — ^ U = u + e(p{t,x,u),e « 1, (8) 

one can check that the invariance analysis of eq.(7) singles out the special value m = 2 
as the only choice having nontrivial Lie point symmetries of the form 

T = a, ^ = b, (p = c{t,x){l-uf, (9) 
where a, b are arbitrary constants and c{t, x) is any solution of the linear heat equation 

Ct-Cxx-c = 0. (10) 
The corresponding symmetry algebra is of the form 

[X,,X2] = 0, = Xe„ [X2,X,] = X,^, (11) 

where 

Xi=dt, X2 = 5xand X^ = c{t,x){l - ufd^, (12) 
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and is of infinite dimensional in nature. However, for all other values of m (to ^ 2) 
one obtains only the trivial translation symmetries 

T = a, ^ = b, (l) = 0. (13) 

Solving now the characteristic equation associated with (9) for to = 2, 

dt dx du 



a b c{t,x){l — uy 

one obtains the similarity variables 



(14) 



z = ax — bt, u=l j:—. r--, (15) 

a + iB{z) + J c{t,x)dt ^ ' 

where w satisfies the similarity reduced ordinary differential equation (ode) of the 
form 

cP'w" + bw + w = 0. (16) 

Eq.(16) is infact the form of the linear heat equation in terms of the wave variable z. 
Further, one can consider the quantity 

1 " 



w{z) + / c{t, x)dt 



X, (17) 



where x satisfies eq.(5), so that the similarity variable u in (15) is nothing but 
the linearizing transformation (6). It has now been given a group theoretical 
interpretation. 

Further, since the general solution of (16) is 



w = he^^^ + he^^^, mi,2 = -^^^^^^, (18) 

where /i and I2 are integration constants, we can write the solution to the original 
pde as 

1 — , 6^- 40^ = 0; 

a + eP(aa^ " bt) (/^ + j^^ax - bt)) + J c{t, x)dt (19) 



u = < 



1 



a + eP("-^ ~ (/i cos q{ax - bt) + h sin q{ax -bt))+J c{t, x)dt ' 

62 - 4a2 < 



with p = —b/2a^, q = \/ ^0? — b'^/2a'^. In the special case c{t,x) = one obtains all 
the interesting travelling wave solutions and stationary structures discussed in [7] . 

For all other values of 2), one obtains from the symmetries (13) the similarity 
variables z = ax — bt and u = w{z) leading to an ode of the form 

a^ww" — ma^w ^ + bwHj' — (1 — w)w^ = (20) 

with w = 1 — w. The above eq.(20) is in general of non-Painleve type when a, 6 ^ 
except for TO = and b = 5/\/6 under proper rescaling [9]. For the other choices 
one can solve for the case 6 = in terms of elliptic function solutions including the 
limiting case of solitary pulse which are all of static form. 

Now we extend our above analysis to the (2+1) dimensional case, 

Ut - Uxx - Uyy - ^ {ul + Uy) - u + u'^ = 0. (21) 
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Under the infinitesimal transformation 

X — * X = X + e£,{t,x,u), y — > Y = y -\- eri{t,x,u), t 



T = t + £T{t, X, u), 



u — >U = u + e(l){t,x,u),e «1, (22) 

the invariance analysis separates out the special value of m = 2 for which eq.(21) 
possesses the following Lie point symmetries; 

r = a, ^ = 632/ + 64, ?7 = -63X + (i4, rp ^ c(t, x,y){l - uf (23) 

where 63, 64 and d4 are arbitrary constants and c{t,x,y) is the solution of the two 
dimensional linear heat equation 

ct - Cxx - Cyy - c = 0. (24) 

But for all other choices of m (7^ 2) we get 

T = a, ^ = 631/ + 64, r] = -b3X + d4, (j) = 0. (25) 

As earlier, for m = 2, the similarity variables arc fonnd by solving the 
characteristic equation associated with the symmetries (23). They are 



zi = — (x^ + t/^) + biy - d4X, Z2 = -t- 



ba 



sm 




w{zi,Z2) + J c{t,x,y)dt 
where w satisfies the similarity reduced (1+1) dimensional pde 



+ "^baw^^ + {2b3Zi + bl + dl)w^^;, 



263Z1 



dl 



+ w ■ 



0. 



(26) 



(27) 



The above equation is nothing but the linear heat equation in terms of the variables 
zi and Z2- Here too as in cq.(15) one can consider the similarity form (26) in two 
dimensions to get the linear heat equation 



Xt Xxx Xyy X — 0) 



1 

x = - 

a 



w{zi,Z2) + / c{t,x,y)dt 



(28) 



so that the transformation for the variable u can again be interpreted as the linearizing 
transformation as in the 1-dimensional case from a group theoretical point of view. 

Carrying out again a Lie symmetry analysis on eq.(27), one obtains the similarity 
variables 

f ClZ2\ 

/(C) + - / C2{zi,Z2)e\ ^3 

C3 J 



\ C3 



-Z2 



dzo 



21=2632:1 + 64 + ^4, ^2=22, 63,^4^0, 

where / satisfies the linear second order ode of the form 

C'/" + C/' + (^ + SC)/ = 0, A = (aci/263C3)2 B = {1 

with prime denoting differentiation w.r.t. C- 

As the general solution to eq.(30) can be expressed in terms of cylindrical functions 
of the form 



ci/c3)/46^ 



(29) 
(30) 



f = hZi{2y/B0+l2Z2{2^B0 



(31) 



Singularity Structure, Symmetries and Integrability 



where Zi's, i = 1,2, are the two linearly independent cylindrical functions and Ii and 
I2 are arbitrary constants, the invariant solution to the (2+l)-diinensional pde (21) is 
written as 



u = 1 — a 



a + e 



ci _ 

C3 



Ii Zi (2 /b^) + /2 (2 /b^) 



C3 J 



ci _ 

Z2 



C2{zi,Z2)e\^3 y dz2\ + / c{t,x,y)dt 



(32) 



Next, proceeding to the special case, 63 = and ^4 = in eq.(26), and carrying out 
an analysis as above, we get the solution to the original pde (21) as 



1 — a< a + exp 



u = { 



-k I —X — t 

M 



h cos(\/ kiC5b4,y) + h sin(v fci 05642/) 



+ j '-^^e^^-dz. 



C5 



1 — a< a + exp 



-k I —X — t 
M 



c{t, X, y)dt 



ki < 0, 



1 — a< a + exp 



, fci>0 



(33) 



+ 



X — t 



hc^b^y + h 



C3{Z1,Z2) kzo. 



where the parameter ki = -J-^ 

blcl 



C5 



k-(— 
^ U4 



dZ2 



c{t, X, y)dt 



fci =0 



with k = —C2/C5 and zi = b^y, Z2 



—X — t. Here C2, C4, C5 are arbitrary constants of integration. 
64 

Eqs.(32) and (33) contain a large class of interesting solutions of various types 
including travelling wave solutions and static patterns. Particularly, we can easily 
check that all the five classes of travelling wave solutions discussed in [7] can be 
derived from eq.(33) for particular choices of the constants involved along with the 
specific choice that the functions 03(^1, ^;2) = and c{t,x,y) = 0. Specifically, the 
simplest travelling wave solution 

^ -, ki>0 (34) 



u = l~ 



1 + Aexp 



—k ( —X — t]± vkic^b^y 



can be constructed by assuming either Ji = or /2 = 0. However, if we choose 
h = I2 (7^ 0) in eq.(33), we obtain the V-wave solution 

1 
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On the other hand for I2 = and /ci < in eq. (33), we obtain the wave front 
oscillating in space: 

u = l F 7 — r-i . (36) 



1 + Aexp 



64 



X — t 



The case fci = 0, /i = describes a travelUng plane wave. But when ^ and 
/2 = we get an inhomogeneous solution 

u = l — 7 r-r (37) 



1 + A\y\ exp 



which is nothing but the separatrix solution. Finally one more choice exists for positive 
fci and h = —l2- 

^ • (38) 



u = 1 



1 + Aexp 



-k I —X — t 
M 



|sinh(v^C564y)| 



This is the Y-wave solution. In each of the above solutions A is a positive constant. 
Several static structures can also be obtained as limiting cases of the above solutions 
(32) and (33). 

Proceeding in a similar fashion for m other than 2 in eq.(21), we obtain the 
similarity reduced variables 



zi = ^(a;^ + 2/^) + b^y - d^x, z2 = -t- sin ^ 
/ 03 

U = w{zi,Z2), (&3, 1^4 7^0). 

along with the reduced pde of the form 

+ "^bsWzi + {2b3Zi + bl + dl)wzizi + 



di — 63X 



^dl + 2b3Zi + bl^ 



(39) 



(2632:1 + bl + dl)wl 



9 9 



2biZi +bl + di 
+ w — vP' = 0. 



w 



' 2bsz, + bl + dj 
On carrying out a similarity reduction. eq.(40) reduces to an ode 



+ f-f = 0, ' = d/dC 



(40) 



(41) 



where 



C = zi, w = f{C), zi=2b3Zi + bl + dl, 



(42) 

giving rise to static structures in [x, y) variables. This equation is found to be 
nonintegrable, in general. However, if we consider for the special case 63 = 0, di = 0, 

the similarity variables become zi = 64J/, Z2 = —x — t, u = w{zi, Z2), eq.(40) on one 

04 

more reduction gives the ode 



Df" 



/''-ci/' + /(l-/) = 0. 



1-/ 
,2 



D=(^^^c^^+bwy '=d/dc. 



(43) 
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with C, = —ci{f^x — t)+ C2hiy and w = /(C), giving rise to plane wave solutions. This 
equation is of non-Painleve type, except for the choice m = in which case eq.(43) 
reduces to the form (20) with m = 0. 

On the other hand with the choice 63 = alone, the similarity variables become 
zi = diX — b^y, Z2 = ax — b^t and u = 'w{zi,Z2). Correspondingly eq.(40) on a further 
reduction reduces to an ode 

+ Bfi - ^f[' - /i + /2 = 0, ' = d/dC (44) 

with /i = 1 - /, A = a'^{cldl + clbf), B = -^464(01 + C2) and C = ac^idiX - 
b^y) — d4{ci +C2)(ax — b/^t), w = /(C)- However, the system is found to possess elliptic 
function solutions including the limiting case of the solitary pulse for particular choices 
of the constants involved. More details of these results will be published elsewhere. 

Throughout our above analysis we have made use of the computer program 
MUMATH [10] to determine the symmetries. 

To conclude, the generalized Fisher type equation with modified diffusion in (1+1) 
and (2+1) dimensions has been foimd to be intcgrable only for the special case m = 2 
via both the singularity structure and symmetry analysis. Moreover we have also 
shown that a Backlund transformation gives rise to the linearizing transformation to 
the heat equation for the integrable case. Besides, a Lie symmetry analysis leads 
in a natural way to the linearizing transformation and physically important classes 
of solutions through similarity variables and similarity reductions, thereby giving a 
group theoretical understanding of the system. 
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